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On the one hand, the explicit Euler scheme fails to converge 
strongly to the exact solution of a stochastic differential equation 
(SDE) with a superlinearly growing and globally one-sided Lipschitz 
continuous drift coefficient. On the other hand, the implicit Euler 
■ scheme is known to converge strongly to the exact solution of such an 

SDE. Implementations of the implicit Euler scheme, however, require 
additional computational effort. In this article we therefore propose 
d ■ an explicit and easily implementable numerical method for such an 

SDE and show that this method converges strongly with the standard 
order one-half to the exact solution of the SDE. Simulations reveal 
that this explicit strongly convergent numerical scheme is consider- 
ed 1 ably faster than the implicit Euler scheme. 
>' 

. 1. Introduction and main result. The explicit Euler scheme (see, e.g., 

Kloeden and Platen [20], Maruyama [26] and Milstein [27]) is most com- 
monly used for approximating stochastic differential equations (SDEs) with 
globally Lipschitz continuous coefficients. Unfortunately, the explicit Euler 
scheme does not converge in the strong mean square sense to the exact solu- 
tion of an SDE with a superlinearly growing and globally one-sided Lipschitz 
continuous drift coefficient. Even worse, Theorem 1 in [16] shows for such 
an SDE that the absolute moments of the explicit Euler approximations at 
K> ! a finite time point T E (0, oo) diverge to infinity. The implicit Euler scheme 
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is better than the explicit Euler scheme in that it converges strongly to the 
exact solution of such an SDE (see Higham, Mao and Stuart [12]). However, 
additional computational effort is required for its implementation. Therefore, 
we wish to identify explicit numerical methods which are strongly conver- 
gent even for SDEs with superlinearly growing coefficients. In this article we 
propose a "tamed" version of the explicit Euler scheme in which the drift 
term is modified such that it is uniformly bounded. Being almost identical 
to the explicit Euler method, this version is explicit and easy to implement. 
Now the benefit of this "tamed" Euler scheme is that it converges strongly 
to the exact solution in case of SDEs with superlinearly growing coefficients. 
More precisely, the main result of this article shows that this "tamed" Euler 
scheme converges strongly with the standard convergence order ^ to the ex- 
act solution of the SDE if the drift coefficient function is globally one-sided 
Lipschitz continuous and has an at most polynomially growing derivative. 
The diffusion coefficient is assumed to be globally Lipschitz continuous here. 
Simulations confirm our theoretical results. 

Throughout the whole article we assume that the following setting is 
fulfilled. Let T S (0, oo) be a fixed real number, let (Q,J-,W) be a prob- 
ability space with normal filtration (J~t)t£[o,T] > let d, m £ N := {1, 2, . . .}, 

let W = {W^\. . . , W (m) ) : [0,T] x O -)• R m be an m-dimensional standard 
(J r i)i e ro i T]-Brownian motion and let £ : S7 — > M. d be an J 7 o/B(M d )-measurable 
mapping with E[||£|| p ] < oo for all p G [l,oo). Here and below we use the 

notation \\v\\ := (\vi\ 2 + • • • + (ffel 2 ) 1 ^ 2 , (v, w) := V\ ■ W\ H + ■ Wk for all 

v = (ui,.. .,v k ), w = (wt,. . .,w k ) €R k ,k£ N, and \\A\\ := sup veR i^ v \\< 1 \\Av\\ 

for all A e R kxl , k, I 6 N. Moreover, let fi : R d -> R d be a continuously differ- 
entiable and globally one-sided Lipschitz continuous function whose deriva- 
tive grows at most polynomially and let a = (cij)ie{i,2,...,d},je{i,2,...,m} :M d — >• 
R be a globally Lipschitz continuous function. More formally, suppose 
that there is a real number c £ (0, oo) such that ||//(x)|| < c(l + ||x|| c ), 
||cr(a;) — <r(y)\\ < c\\x — y\\ and (x — y,/j,(x) — /i(y)) < c\\x — y\\ 2 for all x,y £ 
Then consider the SDE 

(1) dX t = fjL(X t )dt + a(X t )dW t , X = C, 

for t € [0,T]. The drift coefficient \x is the infinitesimal mean of the process X 
and the diffusion coefficient a is the infinitesimal standard deviation of the 
process X. Under the above assumptions, the SDE (1) is known to have 
a unique strong solution. More formally, there exists an adapted stochastic 
process X : [0, T] x Q — > M. d with continuous sample paths fulfilling 

(2) X t = Z+ [ fi(X s )ds+ [ a(X s )dW s 

Jo Jo 

for all t G [0,T] P-a.s. We refer to Theorem 2 in Alyushina [1], Theorem 1 
in Krylov [21] and Theorem 2.4.1 in Mao [24] for existence and uniqueness 
results for SDEs of the form (1). 
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The goal of this article is to solve the strong approximation problem (see, 
e.g., Kloeden and Platen [20], Section 9.3) of the SDE (1). More precisely, 
our aim is to find a numerical approximation Y : Q — > M. d which satisfies 

(3) (E[||z T -y|| 2 ]) 1/2 <e 

for a given precision e > and which can be implemented with as little com- 
putational effort as possible. At this point let us comment on the importance 
of solving the strong approximation problem (3). A central motivation for 
studying strong approximations in the sense of (3) is Giles' seminal paper [7] 
(see also Heinrich [10]). There he introduces, in comparison to the classical 
Monte Carlo method, a very efficient, somehow accelerated Monte Carlo 
method for approximating moments or other expectations of functionals of 
the SDE solution via numerical schemes that converge strongly (see also 
Creutzig, Dereich, Miiller-Gronbach and Ritter [5] for an detailed compari- 
son of the classical and the new "accelerated" Monte Carlo method) . In view 
of this method, strong approximations of the exact solution of the SDE (1) 
in the sense of (3) yield very efficient approximations of expectations of 
functionals of the SDE solution and this is a central reason for developing 
strongly convergent numerical methods. 

The simplest and most obvious idea to solve the strong approximation 
problem (3) is to apply the explicit Euler scheme to the SDE (1). More 
precisely, the explicit Euler method for the SDE (1) is given by mappings 
: n -»• R d , n G {0, 1, . . . , N}, N G N, which satisfy Y N = £ and 

(4) Y n N +1 = Y» + | ■ //(yf) + a(Y n N )(W (n+1)T/N - W nT/N ) 

for all n G {0, 1, . . . , N — 1} and all N G N. In the literature (see, e.g., The- 
orem 10.2.2 in Kloeden and Platen [20], Theorem 1.1 in Milstein [27] or 
Theorem 3.1 in Yuan and Mao [38]) the convergence results for the explicit 
Euler scheme require the drift coefficient \i of the SDE (1) to be globally 
Lipschitz continuous or to grow at most linearly, which we have not assumed 
in our setting. As it turns out, the assumption of an at most linearly growing 
drift function is essentially necessary. More precisely, in the case d = m = 1 , 
it has recently been shown in [16] that the root mean square distance of the 
exact solution of the SDE (1) and of the explicit Euler approximation (4) 
diverges to infinity 

(5) lim (E[||A>-Y^|| 2 ]) 1/2 =oc, 

if the drift coefficient fi of the SDE (1) grows superlinearly, that is, if there 
are real numbers a,C G (l,oo) such that \fi(x)\ > ^Jr- holds for all \x\ > C. 
Thus the explicit Euler scheme (4) does not solve the strong approximation 
problem (3) of the SDE (1) in general. This is particularly unfortunate as 
SDEs with superlinearly growing coefficients are quite important in appli- 
cations (see, e.g., [2, 8, 17, 18, 23, 30, 33]). We remark that in contrast to 
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strong mean square convergence, pathwise convergence of the explicit Eu- 
ler method (4) to the exact solution of the SDE (1) holds due to Gyongy's 
result [9]. 

Another idea for solving the strong approximation problem (3) is to apply 
the implicit Euler scheme, a.k.a. backward Euler scheme (see Higham, Mao 
and Stuart [12]), to the SDE (1). The implicit Euler scheme for the SDE (1) 

is given by mappings ■ ft — > n € {0, 1, . . . , N}, iV£N, which satisfy 

Y$ = £ and 

(6) Y» +1 = Y% + I • M^+i) + a(Y»)(W {n+l)T/N - W nT/N ) 

for all n £ {0, 1, . . . , JV — 1} and all iV£N. A solution of this implicit equation 
is guaranteed to exist and to be unique for N € N large enough due to the 
globally one-sided Lipschitz continuity of /i. In the same setting as in this 
article, Higham, Mao and Stuart showed in Theorem 5.3 in [12] (see also 
[11, 14, 22, 25, 34-37] and the references therein for more approximation 
results on implicit numerical methods for SDEs of the form (1)) that the 
implicit Euler scheme (6) converges with order ^ to the exact solution of 
the SDE (1) in the root mean square sense, that is, they established the 
existence of a real number C € [0,oo) such that 

(7) (EOiXT-y^H 2 ]) 1 / 2 <c- n- 1 ' 2 

for all JVgN. However, additional computational effort is required in order 
to implement (6) since the zero of a nonlinear equation has to be determined 
in each time step in (6). 

To sum up, the explicit Euler scheme, on the one hand, is explicit and 
easily implemented but does, in general, not converge strongly to the ex- 
act solution of the SDE (1). The implicit Euler scheme, on the other hand, 
converges strongly to the exact solution of the SDE (1) but additional com- 
putational effort is required for its implementation. Therefore, we aim at 
a simple explicit numerical method which converges strongly to the exact 
solution of the SDE (1). 

More formally, the following numerical method for approximating the 
solution of the SDE (1) is proposed here. Let Y£ : ft R d , n £ {0,1, ... , JV}, 
N G N, be given by Y^ = £ and 

(8) Y n N +1 = Y n N + Y^T[N^Wl)\\ + a{Y " ){W ^ T / N ~ WnT ' N) 

for all n £ {0, 1, . . . , N — 1} and all N €N. We refer to the numerical meth 

T 

JV 

_ . II ,,(vN-\\\-\ f m „ c /n 1 j 

N 



od (8) as a tamed Euler scheme. In this method the drift term j£ • n(Yj^) 
is "tamed" by the factor 1/(1 + £ • \\fj,(Y^)\\) for n e {0, 1, . . . , N - 1} and 
JV 6 N in (8). Note that the norm of f -fi(Y^)/(l + % ■ \\fJt(Y^)\\) is bounded 
by 1 for every ne{0,l,...,JV— 1} and every N GN. This prevents the drift 
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term from producing extraordinary large values. Additionally, the Taylor 
■ fj,(x) plus terms ^ wuu ^\jp 



expansion of ^ • fi(x)/(l + ^ • ||/i(x)||) in 1/N for fixed x G R is equal to the 
drift term • ^(x) plus terms of order O(j^j). More formally, we see that 

T 



Y» +1 = Y n N + -- nOtf) + a(Y n N )(W (n+1)T/N - W nT/N ) 
(9) 

T\ 2 f,(Y n N ).y(Y n N )\\ 



Nj l + T/N-MY n N )\\ 

for all n G {0, 1, . . . , N — 1} and all JVgN. Thus the tamed Euler scheme (8) 
coincides with the explicit Euler method (4) up to terms of second order. 
Moreover, note that the tamed Euler scheme (8) can be simulated easily and 
the drift function ji needs to be evaluated only once in each iteration of (8). 
More precisely, having calculated v := jj ■ n(Y^), the drift term in (8) is 
then readily computed as i+\\ v \\ ■ 

In order to formulate our convergence theorem for the tamed Euler meth- 
od (8), we now introduce appropriate time continuous interpolations of the 
time discrete numerical approximations (8). More formally, let Y N : [0, T] x 
tt^R d , iVeN, be a sequence of stochastic processes given by 

(10) Y t =Y n+ YTrjN^m\\ i n )( * " nT/N) 

for all t G ne {0,1, .. . ,N- 1} and all iVeN. Note that Y N : 

[0, T] x U — > M. d is an adapted stochastic process with continuous sample 
paths for every N G N. We are now ready to formulate the main result of 
this article. 

Theorem 1.1 (Main result). Let the setting in this section be fulfilled. 
Then there exists a family C p G [0, oo), p G [l,oo), of real numbers such that 

(11) (e[su P \\X t -Y t N \\ p ]) 1/P KCp-N' 1 / 2 

for all N G N and allp G [1, oo). Here X : [0, T] x 17 — > M. d is the exact solution 
of the SDE (1) and Y N : [0, T] xSl4R d , N G N, are the time continuous 
interpolations (10) of the numerical approximations (8). 

Inequality (11) shows that the time continuous tamed Euler approxima- 
tions (10) converge in the strong L p -sense with the supremum over the time 
interval [0,T] inside the expectation to the exact solution of the SDE (1) 
with the standard convergence order i. For a lower bound of this type of 
convergence, the reader is referred to Theorem 3 in Miiller-Gronbach [29] 
(see also Hofmann, Miiller-Gronbach and Ritter [13]). 
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While the detailed proof of Theorem 1.1 is postponed to Section 3, we now 
outline the central ideas in the proof of Theorem 1.1. The key difficulty in the 
proof of Theorem 1.1 is to establish that the tamed Euler approximations (8) 
satisfy the a priori moment bounds 

(12) sup sup E[||F* \\ p ] < oo 

JVeNne{o,i,...,JV} 

for all p G [l,oo) (see Lemma 3.9 in Section 3 for the precise statement of 
this result). After having verified (12), Theorem 1.1 can, by exploiting (9), at 
least in the case p = 2, be completed analogously to Theorem 4.4 in Higham, 
Mao and Stuart [12] in which strong convergence of the explicit Euler method 
under the assumption of the moment bounds (12) has been established. 
However, note that, in contrast to the tamed Euler approximations (8), the 
explicit Euler approximations (4) fail to satisfy such moment bounds for 
SDEs with superlinearly growing coefficients (see (5) here and Theorem 1 
in [16] for details). It is quite remarkable that changing the explicit Euler 
method by a second order term such as in (9) alters the behavior of the 
numerical method to such an extent. 

Let us now go into details and sketch the central ideas of our proof of 
the moment bounds (12). The key idea here for showing (12) is to introduce 
appropriate stochastic processes that dominate the tamed Euler approxi- 
mations (8) on appropriate subevents. More formally, let :fi — > [0, oo), 
rag {0,1,..., iV}, TVgN, be defined by 

/ n-1 

r^:=(A + ||£||)exp(A+ sup £ 



(13) 



A || AW,* 

ue{0,l,...,n}£ 



k 



Y k a r k) A TJ/iV 



\\y k N \\>-i}\ iii^ r j|'"p^ r i 



for all n G {0, 1, . . . , N} and all N G N where A G [1, oo) and AW* : U -> 
R m are defined through A := (1 + 2c + T+ ||/i(0)|| + ||o-(0)||) 4 and AW* := 
W( n+ i) T/JV - W nT/N for all n G {0, 1, . . . ,N - 1} and all N G N. We will 
refer to , n G {0, 1, . . . , N}, N G N, as dominating stochastic processes. 
Appropriate subevents are Sl%, G {0, 1, . . . ,N}, iV GN, given by 



n 



Lei]: sup D%(u) < N 1 '^, 

fc£{0,l,...,n-l} 



(14) 



sup ||AW*(a;)|| <l} 

fc6{0,l,...,n-l} } 



for all n G {0, 1, . . . ,N} and all N G N. The main step of our proof of the 
moment bounds (12) will be to establish the pathwise inequality 

(15) l n M\\Y n N \\<D% 
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for all n G {0, 1, . . . , N} and all N G N (see Lemma 3.1). The next step is 
then to obtain the moment bounds 



(16) limsupE 

7V->oo 



sup \D%\P 

n£{0,l,...,N} 



< oo 



for all p€ [l,oo) for the dominating stochastic processes (see Lemma 3.5). 
These follow nicely from Doob's submartingale inequality (see, e.g., The- 
orem 11.2 in Klenke [19]), from uniform boundedness of a(Y^)/\\Y k , || on 
{\\Y k N \\ > 1} for all k G {0, 1, . . . ,N - 1}, TV G N and from the fact that 



/ 2\Tm T\ 



2 



(17) E[exp(A||A^r + ( U! A^ v ))]<exp^ n ^ + 

for all v G W m , k G {0, 1, . . . , N - 1} and all N G N with N > ATX (see Lem- 
mas 3.2-3.4 for details). Combining (15) and (16) shows that 

(18) sup sup E[l^||y n 7V |n<oo 
JVeNne{0,l,...,AT} 

for all p G [l,oo). For proving (12), it thus remains to verify that 

(19) sup sup E[1 (si n )0 \\Y* \\ p ] <oo 
NeNne{0,l,...,N} 

for all p G [l,oo). This can be achieved by exploiting that the probability 
of (£^) c decays rapidly to zero as N goes to infinity (see Lemma 3.6 for 

T/N-u(Y N ) 

details) and by using that the norm || 1+ jyjv-||^(y^)|| 1 1 — ^ °^ ^ e drift term 
in (8) is bounded by 1 due to the taming factor 1/(1 + ^ • 11/^(^^)11) for 
n G {0, 1, . . . , N} and iY£N. This is, in fact, the only argument in our proof 
of Theorem 1.1 for which the taming factor in (8) is needed. 

It remains to motivate the pathwise inequality (15). The following es- 
timate might at least give an intuition for the case of globally Lipschitz 
continuous coefficients. Consider the tamed Euler scheme for a geometric 
Brownian motion [SDE (1) with d = m = 1, n(x) = and cr(x) = x for all 
iGR], Using 1 + x < e x for i£l yields that 

llyjV u2 = u Y N + Y N AW N f 



(20) 



\Y n N f + liy^AW^H 2 + 2(r»,Y* r AW») 



\Y n N \\ 2 1 + HAWAII 2 + 2( ^W, A_AW? 



yN y N 
1 n 1 n 

\\y»\\' WW 

< \\Y n N fexvl IIA^II 2 + 2( T&m*W» 



! (rai'rar 

for all n G {0, 1, . . . , N — 1} and all JVsN. Iterating this inequality leads to 
an inequality such as (15). Of course, the case of superlinearly growing drift 
coefficients is more subtle and we refer to Section 3.1 for the detailed proof 
of inequality (15). 
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Having sketched the central ideas for the proof of our main result, we now 
compare Theorem 1.1 with some related results in the literature. A priori 
moment bounds of the form (12) for numerical methods for SDEs with non- 
globally Lipschitz continuous coefficients have been intensively studied in the 
literature in recent years. Whereas the references [12, 14, 25, 34-36] deal with 
a priori moment bounds of the form (12) for implicit methods, references for 
explicit methods are infrequent. In particular, in 2002, Higham, Mao and 
Stuart formulated in [12], page 1060, the following open problem: "In gen- 
eral, it is not clear when such moment bounds can be expected to hold for 
explicit methods with f,g£ C 1 " (drift and diffusion coefficients are denoted 
by / and g in [12] instead of [i and a here). In 2010 an a priori moment bound 
for an explicit method has been proved in Bou-Rabee, Hairer and Vanden- 
Eijnden [4]. More precisely, in the setting of the Langevin equation (see 
Section 2 in [4] for the precise assumptions), Lemma 3.5 in [4] establishes ex- 
ponential moment bounds for a version of the metropolis-adjusted Langevin 
algorithm (MALA) with reflection on the boundaries of certain compact 
sets. In this paper we concentrate on the explicit numerical method (8) in 
the setting of the SDE (1). To be more precise, in Lemma 3.9 in Section 3 
below we show that the explicit method (8) fulfills (12) for all p£ [l,oo). 
Furthermore, we observe that one way of deriving the tamed Euler meth- 
od (8) is to approximate the drift coefficent n(v), v E M. d , in the SDE (1) 
by the modified drift coefficients v \\i , v G f° r small A € (0, oo) and 

then to apply the explicit Euler method to these modified SDEs. Approxi- 
mations of this type have been used in the literature in order to construct 
solutions of nonlinear parabolic unilateral problems (see, e.g., Palmeri [31], 
Section 3). Moreover, in the setting of the Langevin equation, Roberts and 
Tweedie suggested in 1996 a similar approximation step to (8) as a proposal 
for the Metropolis-Hastings method in order to sample from the invariant 
measure of the Langevin SDE (see [32], Subsection 1.4.3, (12)-(13)). The 
resulting Metropolis- adjusted method has been named Metropolis- adjusted 
Langevin truncated algorithm (MALTA). Finally, in 1998 a related class of 
numerical methods has been considered in Milstein, Platen and Schurz [28] 
(see [28], (3.2)-(3.3)) in the case of globally Lipschitz continuous coefficients 
of the SDE. To sum up, in the general setting of the SDE (1), the tamed 
Euler method (8) is, to the best of our knowledge, the first explicit numerical 
method that has been shown to converge strongly to the exact solution of 
the SDE (1). 

We now compare simulations of the implicit Euler scheme (6) and of the 
tamed Euler scheme (8). For this we choose T = d = m = 1, £ = 1, /u(x) = 
— x 5 and <t(x) = x for all i£R. The SDE (1) thus reads as 



(21) 



dX t = -Xfdt + X t dWt, X = l, 
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Fig. 1. Root mean square approximation error QE||Xt — Y% || 2 ) 1//2 of the exact solution 
of the SDE (21) and of the implicit Euler scheme (6) and root mean square approximation 
error (E||Xt — Y$ | 2 ) 1 '' 2 of the exact solution of the SDE (21) and of the tamed Euler 
scheme (8) as function of the runtime when N £ {2 4 ,2 5 , . . . ,2 18 }. 

for t £ [0,1]. Suppose that the strong approximation problem (3) of the 
SDE (21) should be solved with the precision of say three decimals, that is, 
with precision e = in (3). Figure 1 depicts the root mean square approx- 
imation error (E||X T -Y^]] 2 ) 1 / 2 of the exact solution of the SDE (21) and of 
the implicit Euler scheme (6) and the root mean square approximation error 
(E||X T - Y$ || 2 ) 1/2 of the exact solution of the SDE (21) and of the tamed 
Euler scheme (8) as function of the runtime when N € {2 4 , 2 5 , . . . , 2 18 }. The 
zero of the nonlinear equation that has to be determined in each time step 
of the implicit Euler scheme (6) is computed approximatively through the 

function fzero(. . .) in Matlab. It turns out that Y 2 H in the case of the 
implicit Euler scheme (6) and that Y 2 ^ in the case of the tamed Euler 
scheme (8) achieve the precision e = j^qq in (3). Following is our Matlab 

code for simulating the implicit Euler approximation Y^t [see (6)] for the 
SDE (21): 

Y = 1; N = 2~16; 
for n=l:N 

v = Y + Y*randn/sqrt (N) ; 

Y = fzero(@(x)x + x~5/N - v, Y) ; 
end 
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Next we specify our Matlab code for calculating the tamed Euler approx- 
imation Y%lt [see (8)] for the SDE (21): 

Y = 1; N = 2~16; 
for n=l:N 
v = -Y"5/N; 

Y = Y + v/(l+abs(v)) + Y*randn/sqrt (N) ; 
end 

The above Matlab code for calculating the implicit Euler approximation Y^it 
requires, on our computer running at 2.33 GHz, about 35.1 seconds while 
the above Matlab code for calculating the tamed Euler approximation Y 2 16 
requires about 0.0212 seconds to be evaluated on the same computer. Thus, 
on the above computer, the tamed Euler scheme (8) for the SDE (21) is 
more than one thousand times faster than the implicit Euler scheme (6) in 
achieving a precision of three decimals in (3). 

2. Further examples. In this section we present further simulations to 
illustrate the efficiency of the tamed Euler scheme (8). The next example is 
a one-dimensional stochastic Ginzburg-Landau equation with multiplicative 
noise (see, e.g., Kloeden and Platen [20], equation (4.52)). More formally, let 
T = d = m = 1, £ = 1, p>(x) = x — x 3 and <r(x) = x for all i£M. The SDE (1) 
thus reads as 

(22) dX t = (X t - X?) dt + X t dW t , X = 1, 

for t € [0, 1]. We now use the implicit Euler scheme (6) and the tamed Euler 
scheme (8) for approximating the SDE (22) and we assume again that the 
strong approximation problem (3) of the SDE (22) should be solved with 
the precision of three decimals, that is, with precision e = j^qq in (3). For 
implementing the implicit Euler scheme (6) in the case of the SDE (22), 
we observe that the drift coefficient of the SDE (22) is a one-dimensional 
polynomial of degree three. Roots of one-dimensional polynomials of degree 
three are known explicitly thanks to Cardano's method. This explicit knowl- 
edge results in a faster implementation of the implicit Euler scheme (6) than 
using the Matlab function fzero(. . .). More precisely, if p, q G R and if q 2 + 
(p/3) 3 > 0, then the only real- valued root of x 3 +px — 2q = 0, x £ R, is x = 

(D + q) 1 / 3 -(D-q)(W where D = y/q 2 + (p/3) 3 £ [0,oo). Thus the implicit 
Euler scheme (6) for the SDE (22) becomes an explicit scheme and satisfies 

Y» +1 = (^) 2 + (iV-l) 3 /27 + ^) 1 / 3 

(23) 

-(VW 2 + (^-l) 3 /27-0 1/3 
for all n 6 {0, 1, . . . , N - 1} and all N G N where q% : U ->■ R is defined by 
q» := ■ f ■ (1 + W {n+1)T/N - W nT/N ) for every n € {0, 1, . . . , N - 1} and 
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10" 5 10"" 10" 3 10" 2 10" 1 

Runtime in seconds 



Fig. 2. Root mean square approximation error (E||AV — Yjvll 2 ) 1 ^ 2 of the exact solution 
of the SDE (22) and of the implicit Euler scheme (23) and root mean square approximation 
error (E\\X T - Y$ | 2 ) 1/2 of the exact solution of the SDE (22) and of the tamed Euler 
scheme (8) as function of the runtime when N £ {2 4 ,2 5 , . . . ,2 18 }. 

every TV € N. Figure 2 displays the root mean square approximation error 

(E||Xr - ^|| 2 ) 1/2 of the exact solution of the SDE (22) and of the implicit 
Euler scheme (23) and the root mean square approximation error (E||Xr — 
Y N ll 2 ) 1/2 of the exact solution of the SDE (22) and of the tamed Euler 
scheme (8) as function of the runtime when TV S {2 4 , 2 5 , . . . , 2 18 }. Comparing 
Figure 1 with Figure 2 confirms that using explicit knowledge of the roots 
of the involved implicit equation results in a much faster implementation 
of the implicit Euler scheme. The tamed Euler scheme, however, is still 
faster as its implementation does not require the arithmetical operations for 

calculating roots. More precisely, it turns out that Y^ll in the case of the 
implicit Euler scheme (23) and that Y?i7 in the case of the tamed Euler 
scheme (8) achieve the desired precision e = in (3). Following is our 

Matlab code for simulating the implicit Euler approximation Y^ 17 for the 
SDE (22) [see (23)]: 

Y = 1; N = 2~17; v = (N-l)~3/27; 
for n=l:N 

q = Y*N*(l+randn/sqrt(N))/2; 

D = sqrt(q"2+v) ; 

Y = (D+qr(l/3) - (D-q)-(l/3); 
end 
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Next we specify our Matlab code for calculating the tamed Euler approx- 
imation Y 2 ^ [see (8)] for the SDE (22): 

Y = 1; N = 2~17; 
for n=l:N 

v = (Y-Y~3)/N; 

Y = Y + v/(l+abs(v)) + Y*randn/sqrt (N) ; 
end 

The above Matlab code for calculating the implicit Euler approximation Y 2 ] 
requires, on our computer running at 2.33 GHz, about 0.0836 seconds while 
the above Matlab code for calculating the tamed Euler approximation Y^ 17 
requires about 0.0467 seconds to be evaluated on the same computer. Thus 
the tamed Euler scheme (8) is on our computer even in the case of the 
SDE (22), where the implicit Euler scheme can be computed explicitly, 
almost two times faster than the implicit Euler scheme (23) in achieving 
a precision of three decimals in (3). 

Our last example is a multi-dimensional Langevin equation. More pre- 
cisely, we consider the motion of a Brownian particle of unit mass in the 
d-dimensional potential z||x|| 4 — ^||x|| 2 , x G where d E N. The corre- 
sponding force on the particle is then x — ||x|| 2 • x, £ G More formally, let 
T = 1, m = d E N, £ = (0, 0, . . . , 0), fi(x) = x — \\x\\ 2 ■ x and let a(x) = I be the 
identity matrix for all x E Mr. Thus the SDE (1) reduces to the Langevin 
equation 

(24) dX t = (X t -\\X t \\ 2 -X t )dt + dW u X = 0, 

for t E [0, 1] . Here W : [0, 1] x $7 — > M. d is a <i-dimensional standard Brown- 
ian motion. Figure 3 displays the root mean square approximation error 

(E\\X T - Y^f) 1 / 2 of the exact solution of the SDE (24) with d = 10 and 
of the implicit Euler scheme (6) and the root mean square approxima- 
tion error (E||X T - *jv || 2 ) 1/2 of the exact solution of the SDE (24) with 
d = 10 and of the tamed Euler scheme (8) as function of the runtime when 
N E {2 4 , 2 5 , . . . , 2 18 }. We see that both numerical approximations of the 
SDE (24) apparantly converge with rate 1. This is presumably due to the 
additive noise in (24). Note that we used the Matlab function fsolve(. . .) in 
our implementation of the implicit Euler scheme (6) for the SDE (24) as the 
Matlab function fzero(. . .) used for the numerical simulations in Figure 1 
is restricted to one dimension {d = 1). 

Some applications involve high-dimensional SDEs (see, e.g., Beskos and 
Stuart [3], Section 2.1). Then the above implementation of the implicit Euler 
method has an additional disadvantage. The Matlab command fsolve(. . .) 
uses (by default) the "trust-region-dogleg" algorithm which calculates Ja- 
cobian matrices. Thus the computational effort increases quadratically with 
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Fig. 3. Root mean square approximation error {E\\Xt 



2\l/2 



of the exact solution 



of the SDE (24) with d = 10 and of the implicit Euler scheme (6) and root mean square 
approximation error (E||Xt — Yff || 2 ) 1//2 of the exact solution of the SDE (24) with d = 10 
and of the tamed Euler scheme (8) as function of the runtime when N £ {2 4 , 2 5 , . . . , 2 18 }. 



the dimension d £ N for every fixed N £ N. To visualize this we have plot- 
ted the runtime of the calculation of the implicit Euler approximation i^J^s 
of the SDE (24) as a function of the dimension d G {10, 20, 30, . . . , 400}. 
Figure 4 suggests a quadratic dependence of the runtime of the implicit 
Euler method on the dimension in case of the SDE (24). In contrast, the 
tamed Euler method is linear in the dimension (except that the evaluation 
of the coefficient functions might increase quadratically in the dimension). 
Figure 5 suggests a linear dependence of the runtime of the tamed Euler 
method on the dimension d £ {10, 20, 30, ... , 400} in case of the SDE (24). 

3. Proof of Theorem 1.1. In order to simplify the notation introduced 
in Section 1, we introduce the mappings a 1 ^ : Q, — > M defined by 

(25) a n .= 1 {|| ^ || > 1} ^_,__AW' B J 

for all n € {0, 1, . . . , N — 1} and all iVsN. Using this notation, the domi- 
nating stochastic processes n £ {0, 1, . . . , N}, N G N, [see (13)] simplify 
to 

D% = (A + ||f ||) • exp (A + sup f>|| AW k N f + a$]\ 
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Fig. 4. Runtime of the calculation of the implicit Euler approximation Y\%% of the 
SDE (24) as a function of the dimension d € {10, 20, 30, 400}. 
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Fig. 5. Runtime of the calculation of the tamed Euler approximation Y^i of the 
SDE (24) as a function of the dimension d £ {10, 20, 30, 400}. 

for all n G {0, 1, . . . , N} and all N € N. Moreover, we denote by e% = (1, 0, . . . , 
0) G R m , e 2 = (0, 1, 0, . . . , 0) G R m , . . . , e m = (0, . . . , 0, 1) G R m the unit vec- 
tors in the M m . Additionally, we use the mappings <Tj : R d — >M. d , i £ {1,2, . . . , 
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m}, given by o~i(x) = (a\^(x), . . . , a^i(x)) = <r{x)ei for all x G M. d and all 
i € {1,2, ... , m}. The SDE (2) can thus be written as 



x t =i+ ii{x s )d S +y] / ai(x s )dwp 

Jo i=l Jo 



(26) 



for all t G [0,T] P-a.s. Our proof of Theorem 1.1 relies on the following 
lemmas. 

Lemma 3.1 (Dominator lemma). Let Y^-.Q^ R d , D% : O -)■ [0, oo) and 
£1% G J 7 /or n G {0, 1, . . . , N} and N G N 6e flwen 67/ fi5j and f4 ). Then 
we have that 



(27) ln^ll^ll<A. 
for all n G {0, 1, . . . , N} and all N G N. 

Lemma 3.2. Let n G N and let Z M. n be an n- dimensional standard 
normal random variable. Then we have that 

(28) E[exp(a||Z|| 2 )] = (1 - 2a)" n/2 < e 2an 
for all a G [0, \]. 



N 



Lemma 3.3. We have that 



(29) 



for all p G [1, 00) 



sup E 

N>AXpT 



N-l 



cxp 



fc=0 



AM 



ATii2 



< OO 



Lemma 3.4. Let :Sl 
by (25). Then we have that 



for n G {0, 1, . . . , N} and N G N be given 



(30) 

for all pG [1, 00 



sup sup E 

«6{-l,l}JVeN 



sup exp pz > a k 



ne{0,l,...,N} 



n.-l 



A' 



fc=0 



< OO 



Lemma 3.5 (Uniformly bounded moments of the dominating stochastic 
processes). Let : — > [0, oo) for n € {0, 1, ... , iV} and TV G N be given 
by (13). Then we have that 



(31) 

for all p G [1, oo). 



sup E 



sup \D, 



Nip 



N<=N, L n£{0,l,...,N} 
N>8XpT 



< oo 
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Lemma 3.6 (Estimation of the probability of the complement of fijy for 
N G N). Let Q$ G J for N G N be given by (14). Then we have that 



(32) 

for all p G [1, 00). 



sup(iV p -P[(^) c ]) <oo 

JVeN 



Lemma 3.7 (Time continuous Burkholder-Davis-Gundy type inequality). 

Let k G N and let Z: [0, T] x O — > M fcxm &e a predictable stochastic process 
satisfying P[J Q ||^ s || 2 ds < 00] 

(33) 

for all t G [0, T] and all p G [2, 00) . 



1. Then we obtain that 

t m 











sup 


r z u dw u 




<p([ 


se[o,t] 


Jo 




LP(U;R) \J0 



1/2 



£11* 

i=l 



s e i\\Lv(a,] 



ds 



ikxm 



Lemma 3.8 (Time discrete Burkholder-Davis-Gundy type inequality). 
Let k G N and let Zf-.Sl^ R kxm , I G {0, 1, . . . , N — 1}, N G N, be a family 

Fit/n I B(R kxm ) -measurable for all 

1/2 



is 



of mappings such that Z^ : U 

I G {0, 1, . . . , N - 1} and all N G N. Then we obtain that 

/ n—l m 



(34) 



sup 

je{0,l,-,n} 



/2 Z l* AW i 
1=0 



<p ££ n*. 



LP(Q;R) 



J=0 j=l 



T 



for all n G {0,1,..., N} ; N G N and a// p G [2, 00). 



Lemma 3.9 (Uniformly bounded moments of the tamed Euler approxi- 
mations). Let : Q ->■ ]R d /or n G {0, 1, . . . ,N} and N G N be given by (8). 
Then we have that 



(35) 

for all p G [1, 00) 



sup sup ElUE^f] <oo 
iVeNne{0,l,...,jV} 



Lemma 3.10. Let Y* 1 : Q -> M d 
oy T/ien roe /i<roe £/ia£ 

sup sup E[|| M (y n JV )|H<oo 

JVeNne{0,l,...,7V} 

/or a// p G [1, 00). 



for n G {0, 1, . . . , N} and N G N be given 



sup sup E[||o-(E n 7V )|H<oo 
Nmn£{0,l,-,N} 



The proofs of Lemmas 3.1-3.10 can be found in Sections 3.1-3.10. Us- 
ing Lemmas 3.9 and 3.10, the proof of Theorem 1.1 is then completed in 
Section 3.11. 
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3.1. Proof of Lemma 3.1. First of all, note that \\AW^ || < 1 on 
for all n G {0, 1, . . . , N — 1} and all JV£N. The global Lipschitz continuity 
of a and the polynomial growth bound on // (see Section 1) therefore imply 
that 

Kill < \\Y n N \\ + §11/^)11 + \\o-(y»)\\\\aw»\\ 

<l + T\HY n N )-v(0)\\+T\M0)\\ 

+ \\ (r (Y»)-a(0)\\ + \\a(0)\\ 
< l + 2Tc + T||/i(0)|| +c+ \\a(0)\\ < A 
on FL% +1 n{we(l: II^HH < 1} for all n G {0, 1, . . . ,N - 1} and all JVgN. 

2 l2 

Moreover, the Cauchy-Schwarz inequality and the estimate a-b < ^- + 
for all o,kR show that 

2 



l J n+ll 



l + (T/iV) HMWII 



II 



11 » 11 + (i + (r/iv)|| M (y r f)||)2 + 11(71 » ' 11 

+ I rrn/ AT \U. /vAfMI + 2 ( r n .^Vn ) AM n ) 



l + (T/N)\\rtV»)\\ 
2((T/N)v(Y n N ),a(Y n N )AW 1 



l + (T/N)\\»(Y n *)\\ 

(37) < ||y^|| 2 + ^IKOll 2 + Ikttf )fl|Awf || 2 

(2T/jV)(y„ jV , / ,(y ? f)) „ , YN , Aw m 

+ ^ll/*0^)lllk(i^)IIIIAw^|| 

< ll^ll 2 + ^IM^)f + 2||a(^)|| 2 ||AW f ri| 2 

(2T/iv)(y n jV , / ,(y ? f)) „ _^ AI ^ 
+ 1 + ^)11^)11 +2(K " )A ^ } 

on Q for all n G {0, 1, . . . , N - 1} and all N G N. Additionally, the global 
Lipschitz continuity of a (see Section 1) implies that 

\W(x)\\ 2 < (\\a(x) -a(0)|| + |K0)||) 2 

(38) < (c||x|| + ^(O)!!) 2 < (c+ ||(j(0)||) 2 ||x|| 2 
< Allxll 2 
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for all x G M d with ||x|| > 1 and the global one-sided Lipschitz continuity 
of fi (see Section 1) gives that 

(x, fi(x)) = (x, fi(x) - fi(0)) + (x, fi(0)) 

(39) <c||x|| 2 + ||x|||| M (0)||<(c+|| M (0)||)||x|| 2 

< v^Hxll 2 

for all with ||x|| > 1. Furthermore, the polynomial growth bound 

on // (see Section 1) yields that 

H^II 2 <(H^)-M0)ll + H0)ll) 2 

\2 



(40) < (c(l + \\x\\ c )\\x\\ + |K0)||r < (2c\\x\H c+1) + 

< (2c+ ||/i(0)||) 2 ||x|| 2(c+1) < Ny/\\\x\\ 2 

for all x G M. d with 1 < ||x|| < TV 1 /^) and all N G N. Combining (37)-(40) 
and T 2 + T < y/\ then gives that 

ll^ill 2 < \\Y n N f + ^ll^ll 2 + 2A||r r rf HAWAII 2 

+ ^\\Y n N \\ 2 + 2(Y n N MY n N )AW») 

(41) 

= \\Y n N f (l + 2(T ^ T) ^ + 2A|| A^|| 2 + 2a N n 

<rn 7V H 2 exp(^+2A||A^|| 2 + 2^ 

on {oj G 0:1 < < N 1 '^} for all n G {0, 1, . . . , N - 1} and all 

NeN. 

Additionally, we use the mappings rf : O — > {— 1, 0, 1, . . . , I — 1}, I G {0, 1, 
...,N}, N GN, given by 

(42) r/» := max({-l} U {n G {0, 1, . . . , I - 1}\\\Y» ( w )|| < 1}) 

for all cj G O, Z G {0, 1, . . . , N} and all N G N. 

With the estimates (36) and (41) at hand, we now establish (27) by induc- 
tion on n G {0, 1, . . . , N} where N G N is fixed. The base case n = is trivial. 
Now let I G {0, 1,. . . ,N — 1} be arbitrary and assume that inequality (27) 
holds for all n G {0, 1, . . . , I}. We then show inequality (27) for n = I + 1. 
More formally, we now establish that 

(43) < A+iM 

for all oj G fi^L. To this end let to G QfLi be arbitrary. From the induction 
hypothesis and from u G Vtf +1 C it follows that ||Y^(w)|| < D%(u) < 
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N i/(2c) for all n g jo, 1, . . . , Z}. By definition (42) we therefore obtain that 1 < 
HY^HH KN 1 '^ foraUne{7^ 1 (w) + l,7^ 1 (w) + 2,...,Z}. Estimate (41) 
thus gives that 

(44) \\Y n N +1 (u)\\ < ||Y*(«)|| -ex P (^ + A||A^|| 2 + a^)) 

for all n £ {tK^u) + 1, tK^lj) +2,..., I}. Iterating (44) hence yields that 



<---<\\Y N 



< \\y3 



% Auj)+1 ( W ) || • exp ( ]T (4 + A H AW n(") II 2 + «n M) ) 

& M+1 (uO||-exp(A+ sup £(A||AW^( W )|| 2 + a?M)Y 



Estimate (36) therefore shows that 



(45) < (A + ||£( w )||) • exp (x + sup ]T(A|| AW? (co)\\ 2 + a% (V 

\ ue{o,i,...,i+i} n=u j 

This finishes the induction step and the proof. 

3.2. Proof of Lemma 3.2. We establish (28) in the case n = l. The gen- 
eral case then follows from independence. For the case n = 1 note that 

(46) < e 2 *' 

(l-x) 

for all x £ [0, Inequality (46) then shows that 

E[exp(aZ 2 )]= / e ax " -^e^ 2 dx = [ -L e - 2 /(2(l-2a)-i) dx 

1 <e 2a 



VI -2a 

for all a 6 [0, j]. This completes the proof of Lemma 3.2. 
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3.3. Proof of Lemma 3.3. The independence of the random variables 
AWjf , k G {0, 1, . . . ,N — 1}, and Lemma 3.2 show that 



E 



/ JV-1 
exp I Xp ^ 
V fc=o 



\AW k N f 



N-l 



I E[exp(Ap||A^f || 2 )] 



k=0 



(47) 



= (E[exp(Xp\\W T/N \\ 2 )]) N 
< exp(2ApTm) < oo 



for all JVsNfl [4XpT, oo) and all p G [l,oo). This completes the proof of 
Lemma 3.3. 



3.4. Proof of Lemma 3.4- First of all, note that the time discrete stochas- 
tic process z]Tfe=o n G {0, 1, . . . ,N}, is an (J r nT / A r)ne{o,i,.--,A r }" mar ^ n S a l e 
for every z £ { — 1, 1} and every N G N. In particular, we therefore obtain 



that the time discrete stochastic process exp(z^^~ a k ), n G {0, 1, . . . , iV}, 
is a positive (J 7 n T// v r) ne {o,i,...,7V}-submartingale for every z G { — 1, 1} and ev- 
ery iV G N. Doob's maximal inequality (see, e.g., Klenke [19], Theorem 11.2) 
hence shows that 



(48) 



n-l 



sup exp z 

ne{0,l,...,7V} V k=0 



< 



p 



(p-1) 



exp z 



N-l 

£ 

fc=0 



Lp(n-,R) 



for all iV G N, p G (l,oo) and all z G { — 1, 1}- Moreover, we have that 
E 



% o~(x) . 

p z H\\x\\>i}{\UJ\> j^f AW > 



N 

p 2 T 
N 



Ml>i}" 



(a(x))*x\\ 2 < p 2 T 



N 



-1 



{Nl>i}" 



* l|2 



< 



p 2 T(c+ \\a(0)\\Y' 
N 



for all x G M d , fe G {0, 1, . . . , N - 1}, iV G N, p G [l,oo) and all z G {-1,1}. 
Lemma 5.7 in [15] therefore gives that 

_E_ ^a^XM fp 2 T(c+ ||a(0)||) 2 



(49) E 



exp I pzl { || x ||>i } 



< exp 



iV 



for all x G M d , A; G {0, 1, . . . ,N - 1}, iV G N, p G [l,oo) and all zG{-l,l}. 
Estimate (49), in particular, shows that 



E[exp(pzaf )|J r fcT/A r] < exp 



p 2 T{c+ ||<r(0)||)' 
N 
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for all k G {0, 1, . . . ,N - 1}, N G N, pG [l,oo) and all 2 G {-1, 1}. Hence, we 
obtain that 



E 



(50) 



/ JV-l 



exp 



k=0 



E 



< E 



7V-2 



exp I pz ^ of J ■ E[exp(p^ajv_i)|^"(jv-i)T/Jv] 



fc=0 
7V-2 



exp 



fc=0 



• exp 



/ p 2 T(c+ ||a(0)||)- 
V N 



<---<exp(p 2 T(c+\\a(0)\\) 2 ) 

for all iV G N, p G [1, oo) and all z G { — 1, 1}- Combining (48) and (50) then 
gives that 

/ n-l 



sup sup 



exp I z a 
V k=0 



L*(fi;] 



ne{0,l,...,JV} 

< 2exp(p 2 T(c+ |K0)||) 2 ) < oo 
for all p G [2, oo) and this completes the proof of Lemma 3.4. 



3.5. Proof of Lemma 3.5. First of all, Holder's inequality shows that 



sup 

iVeN, 

N>8XpT 



sup D„ 

ne{0,l,...,7V} 



(51) 



<e A (A+ ||£||i4 P(nKd) ) sup 

N>8XpT 



N-l 



exp a£||AW£ 



Ni\2 



k=0 



L 2 p(U;) 



x sup 
\ JVeN 



n-l 



sup exp sup y, 

ne{0,l,...,JV} \we{0,l,...,n} k=u 



L a p(Q;\ 



for all p£ [l,oo). Moreover, again Holder's inequality implies that 



n-l 



sup exp sup a k 



ne{0,l,...,N} 



(52) 



< 



ue{o,i....,»} fc=u 

n-l 

sup exp ( ^ a% 



A? 



ne{0,l,...,JV} 



u-1 



sup exp 

ue{o,i,...,iV} 



fc=0 



L 8 P(Q;I 
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for all N G N and all p G [1, oo). Putting (52) into (51) then gives that 



sup 

JVeN, 

N>S\pT 



sup D% 



ne{0,l,...,N} 



(53) 



<e A (A+||e||L4p(n;K)) 



sup 

iVeN, 

"iV>8ApT 



7V-1 



exp [\J2\\AW£ 



N\\2 



sup sup 

ze{-i,i}iV6N 



k=0 
n-1 



sup exp [ z 

n£{0,l,...,N} 



fc=0 



\ 2 



for all p G [l,oo). Combining Lemmas 3.3, 3.4 and (53) finally completes the 
proof of Lemma 3.5. 

3.6. Proof of Lemma 3.6. The subadditivity of the probability mea- 
sure P[-] and Markov's inequality yield that 



P[(fi 



n) J 



(54) 



< P 

< E 

< E 



sup D^yN 1 '^ 
ne{0,l,...,N-l} 

sup \D*\ q 

L nG{0,l,...,iV-l} 



sup |D r 

ne{0,l,...,N} 



Nw 



+ N-¥[\\W T/N \\ > 1] 

.N-i/W +N-F[\\W T \\ >v / iV] 
Ar-9/(2c) +E [||p FT ||9]. iV (i-g/2) 



for all N G N and all (/ G [l,oo). Combining Lemma 3.5 and (54) then 
shows (32). This completes the proof of Lemma 3.6. 

3.7. Proof of Lemma 3.7. Lemma 3.7 is an immediate consequence of 
Doob's maximal inequality, of the Burkholder-D avis-Gundy type inequality 
in Lemma 7.7 of Da Prato and Zabczyk [6] and of the triangle inequality. 
For completeness we now present the proof of Lemma 3.7. 

Proof of Lemma 3.7. Doob's maximal inequality (see, e.g., Da Prato 
and Zabczyk [6], Theorem 3.8), Da Prato and Zabczyk [6], Lemma 7.7, and 
the triangle inequality give 
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= \P\ P 



v 



2(p-l 

t m 

o 



p/2 / r t m 

En*. 



\ P/2 



liP(n;K fe ) 



i=l 



p/2 



P-ll 2 



ds 



i=i 



for all t G [0,T] and all p G [2, oo). This completes the proof of Lemma 3.7. 

□ 



3.8. Proof of Lemma 3.8. Let Z N : [0,T] x 



R kxm , N G N, be a se- 
quence of stochastic processes defined by := Zf for all s£[S, 
/G{0,l,...,iV-l} and all iV G N. The Burkholder-Davis-Gundy type in- 
equality in Lemma 3.7 then shows 
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for all re G {0, 1, . . . , N}, N G N and all p G [2, oo). This completes the proof 
of Lemma 3.8. 

3.9. Proof of Lemma 3.9. In order to show Lemma 3.9 we first represent 
the numerical approximations (8) in an appropriate way. More formally, we 
have that 

n-1 



y, 
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(T/NMY k N ) 
+ (T/N)MY")\\ 



n-1 



N 
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A 1 
fe 



fc=0 



P-a.s. for all re G {0, 1, ... , iV} and all iV G N. The Burkholder-Davis-Gundy 
type inequality in Lemma 3.8 then gives that 

rN I 



n llLP(n;]E 
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and the global Lipschitz continuity of a therefore shows that 
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for all re 6 {0, 1, ... , iV}, JVgN and all p G [2, oo). In the next step Gronwall's 
lemma gives that 
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ne{0,l,...,JV} 



n \\LP(fl; 



< V2e? 2mTc2 m\ LP{n . iRd) + py/T^\\a(0)\\+N) 



for all N £ N and all p € [2, oo). Of course, (58) does not prove Lemma 3.9 
due to the JVGNon the right-hand side of (58). However, exploiting (58) in 
an appropriate bootstrap argument will enable us to establish Lemma 3.9. 
More formally, Holder's inequality, estimate (58) and Lemma 3.6 show that 
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x (lieili2 P(n . Rd) + 2px^||a(0)|| + l) 



< oo 

for all p€ [2,oo). Additionally, Lemmas 3.1 and 3.5 give that 
sup sup IIIq^^IIlp^^) 

WeN, n£{0,l,...,N} 
N>8XpT 

(60) 

< sup sup \\D^\\ LP{n . R) <oo 

WGN, n€{0,l,...,N} 
N>SXpT 

for all p G [l,oo). Combining (59) and (60) finally completes the proof of 
Lemma 3.9. 

3.10. Proof of Lemma 3.10. Lemma 3.10 follows from Lemma 3.9. More 
precisely, the estimate \\fi(x)\\ < (2c + ||/i(0)||)(l + ||x||( c+1 )) for all x G R d 
and Lemma 3.9 give that 

SUp SUp \\K Y n)\\LP(nm 
N&Nn£{0,l,...,N} 

(61) <(2c+|| M (0)||)fl + sup sup ll^llSci)^ 



NeNnG{0,l,...,N} 



< OO 



for all p £ [l,oo). Additionally, the inequality ||a"(x)|| < c\\x\\ + ||cr(0)|| for all 
x G M. d and again Lemma 3.9 show that 

SUp SUp 11^(^)11^(0:^) 
NeNne{0,l,...,N} 

(62) <c(sup sup \\Y^\\ Lnn . Rd) ) + ||(7(0)|| 

< OO 

for all p G [l,oo). This completes the proof of Lemma 3.10. 

3.11. Proof of Theorem 1.1. Using Lemmas 3.9 and 3.10 we now estab- 
lish inequality (11). To this end we use the notation 

f f T 2T (N-1)T 1 
Wat :=max|se|0, N ,T\:s<t 

for all t G [0,T] and all N G N. In this notation, equation (10) reads as 
(63) 



j 
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for all s E [0, T] P-a.s. and all N E N. Our goal is then to estimate the 
quantity E[sup sg r 0)t i \\X S — Y S N \\ P ] for t E [0, T] and p E [l,oo). To this end 
note that (2) and (63) imply that 



X s -Y s 
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„(X \ PV V u \n' 

\ 1 + (T/N)\MY^ n )\\ 



da 
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+ J2 / (v^-aiiY^dWU 
i=i ^° 

for all s E [0,T] P-a.s. and all N E N. Ito's formula hence gives that 
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and the inequality (a + 6) 2 < 2a 2 + 26 2 for all o,6gR, the estimate a < |a| 
for all a E M and the Cauchy-Schwarz inequality therefore yield that 



\X S - F s iV Y < ( 2c + 2cm ) / \\ x u - yj II du 

Jo 



+ 2/ \\X u -Y«\\MY u ")-»(Y» }N )\\du 
Jo 



2T 



*' 



+ 2 



i=l ^° 



+ 2c 2 ? n/ ll^-yfjjl 2 ^ 



SDES WITH NONGLOBALLY LIPSCHITZ CONTINUOUS COEFFICIENTS 27 



for all s G [0,T] P-a.s. and all N G N. The inequality a • 6 < % + ^ for all 
a, 6 G R then shows that 



sup \\X S - Y t 



AT || 2 
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<2(c + cWl)/ \\X S -Y S N \\ 2 ds + I \\fi(Y^ ) - fi(Y^ N )\\ z ds 
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P-a.s. for all t G [0,T] and all N G N. The Burkholder-Davis-Gundy type 
inequality in Lemma 3.7 hence yields that 
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for all t G [0, T], iV£ N and all p G [4, oo). Next the Cauchy-Schwarz in- 

2 1.2 

equality the Holder inequality and again the inequality a • b < ^- + for all 
a, b G R imply that 
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< I sup \\X S - Y £ 
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for all t£ [0,T], iVsN and &llp£ [4, oo). Inserting inequality (65) into (64) and 
applying the estimate (a + b) 2 < 2a 2 + 2b 2 for all a, b £ K then yields that 
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for all t € [0,T], iV6N and all p G [4, oo). In the next step Gronwall's lemma 
shows that 
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and hence, the inequality \J a + b + c < y^a + v& + y/c for all a, b, c G [0, oo) 
gives that 



A' I 
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< y/2T e 2T ('P 2 c 2 'in+c+l) I 
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for all N G N and all p G [4, oo). Additionally, the Burkholder-Davis-Gundy 
type inequality in Lemma 3.7 shows that 
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for all iV G N and all p G [2, oo). Lemma 3.10 hence implies that 
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for all p G [l,oo). In particular, we obtain that 
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for all p G [1, oo) due to Lemma 3.9. Moreover, the estimate \\fJ,(x) — fJ,(y)\\ < 
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for all N G N and all p G [l,oo) and inequalities (67) and (68) hence show 
that 



for all p G [l,oo). Combining (66), (67), (70) and Lemma 3.10 finally 
shows (11). This completes the proof of Theorem 1.1. 
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